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Abstract 

Using spatial domain techniques developed by the authors and Myunghyun Oh in 
the context of parabolic conservation laws, we establish under a natural set of spectral 
stability conditions nonlinear asymptotic stability with decay at Gaussian rate of spa- 
tially periodic traveling-waves of systems of reaction diffusion equations. In the case 
that wave-speed is identically zero for all periodic solutions, we recover and slightly 
sharpen a well-known result of Schneider obtained by renormalization/Bloch transform 
techniques; by the same arguments, we are able to treat the open case of nonzero 
wave-speeds to which Schneider's rcnormalization techniques do not appear to apply. 

1 Introduction 

In this paper, we study the nonlinear stability with respect to spatially localized pertur- 
bations of spatially periodic traveling- wave solutions u(x,t) = u(x — ct) of a system of 
reaction diffusion equations of form ut = u xx + f(u), where (x,t) € R x R + , u £ R n , and 
/ : R n — > W 1 is sufficiently smooth: equivalently, spatially periodic standing-wave solutions 
u(x, t) = u(x) of 

(1.1) u t - cu x = u xx + f(u). 

For the Allen-Cahn (variational) case f(u) = dF(u), the traveling-wave ODE becomes 
— cv! = u" + dF(u), hence for c ^ either increases or decreases a Hamiltonian 

lu'l 2 

H:= l -^-+F(u), 
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from which it follows that periodics exist only for speed c = 0. In this zero-speed case 
(not necessarily originating from variational form), Schneider [Sj proved nonlinear stability 
with decay at Gaussian (diffusive) rate by a combination of weighted energy estimates, 
renormalization techniques, and a spectacular nonlinear cancellation estimate, all carried 
out in the Bloch frequency domain. In the process, he showed that behavior is purely 
diffusive, with no associated convection. See also [CQ, IMSUj and references therein. 

However, as described in [DSSSj . there exist many other cases for which there exist 
spatially periodic solutions of varying speed, in which situation one expects asymptotic 
behavior driven by nonzero convection as well as diffusion. In this case, the renormalization 
argument of [S] does not directly apply. 

Meanwhile, motivated strongly by the work of Schneider, we and co-authors have carried 
out by somewhat different methods stability of spatially periodic solutions of systems of 
parabolic conservation laws Ut + h(u) x = u xx , for which convection plays a major role 
\OZ\ IJZj . Natural questions are (i) whether these alternative methods might be used to 
reproduce the original results of Schneider in the zero-speed case, and (ii) whether, more, 
they might be able to treat the nonzero-speed (convective) case left open up to now. 

In this paper, we answer both questions in the affirmative, establishing stability and 
decay at Gaussian rate for the general case, with no condition on the wave-speed; see 
Theorem 14. II in Section^ Moreover, our nonlinear iteration method, based on spatial rather 
than (Bloch) frequency domain, permits an extremely brief and simple proof yielding new 
insight even in the zero-speed case treated previously by Schneider. 

Note: We have been informed by Bjorn Sandstede that a similar result has been ob- 
tained by different means in }SSSUl using a nonlinear decomposition of phase and amplitude 
variables as in [DSS S1. accomodating also nonlocalized perturbations in the phased 

Acknowledgement. Thanks to Blake Temple, for his encouragement and interest in 
the problem. Also, we are greatful to Bjorn Sandstede for describing the setting of the 
nonzero speed problem and for pointing out the reference [SSSUJ. 

2 Existence and Spectral Stability Assumptions 

Any standing wave solution of (jl.ip clearly must be a solution of the ordinary differential 
equation 

(2.1) u xx + cu x - f(u) = 0, 

which is commonly referred to as the profile equation or as the traveling wave ODE corre- 
sponding to the original system. The existence of periodic orbits of f)2. 1 1) is trivial in the 
case c = where the equation is clearly Hamiltonian and hence can be directly integrated by 
quadrature. When the wave-speed is non-zero, however, the existence of periodic orbits is 
more delicate but still straightforward: it can be treated via ODE/imp licit function theorem 

1 In the notation of Thm. 14.11 data {to with ||wo(ai + V'o) — wILintynH* (R) Ij \®"''1 > o\h k (K) sufficiently small. 
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techniques as familiar in the conservation law case. Indeed, writing (|2.ip as a 2n x 2n sys- 
tem with In constraints (periodicity) and two extra parameters (speed and period) yields, 
generically, a two-dimensional solution set. The techniques presented in this paper apply 
regardless of whether c is non-zero and hence we study (|2.ip for a general wave-speed cGK. 

Throughout our analysis, we assume the existence of an X-periodic solution u(x) of 
(|2.ip . It follows that generically one expects the periodic orbits to form a two-parameter 
family of solutions of the profile equation, parameterized by the wave-speed c and a trans- 
lation mode Xq. More precisely, we make the following generic assumptions: 

(HI) / € C K (M) for some K > 2. 

(H2) The set of periodic solutions of (]2.ip in the vicinity of u forms a smooth 2-dimensional 
manifold {u xo,s (x — st + xq)} with xo, s S IR 

We begin our stability analysis by considering the linearization of (jl.ip about the fixed 
periodic standing- wave solution u. Without loss of generality, we assume that u is 1-periodic, 
i.e. that u(x + 1) = u(x) for all iGl. Considering nearby solutions of the form 

u(x) + ev(x,t) + 0(e 2 ), 

where \e\ <C 1 and v(-,t) £ L 2 (M), corresponding to spatially localized perturbations, we 
see that v satisfies the linear equation 

v t = v xx + df{u)v. 

Since this equation is autonomous in time, we may seek separated solutions of the form 

v(x,t) = e^v(x) 
which readily yields the spectral problem 

(2.2) fiv = L[u}v := {d 2 x + cd x + df(u)) v 

considered on the real Hilbert space L 2 (M), where here we are considering the linear operator 
L[u] as having dense domain in H 2 (W). As coefficients of L[u] are 1-periodic, Floquet theory 
implies the I? spectrum is purely continuous and corresponds to the union of the L°° 
eigenvalues corresponding to considering the linearized operator with boundary conditions 
v (x + T) = e tK v(x) for all where k € [— vr,7r] is referred to as the Floquet exponent 

and is uniquely defined mod 2tt. In particular, u £ a(L[u]) if and only if the spatially 
periodic spectral problem (|2.2p admits a bounded eigenfunction of the form 

(2.3) v(x) = e^ x w(x) 
where w(x + 1) = w(x). 
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Substitution of the Ansatz (|2.3p into (|2.2p motivates the use of the Fourier-Bloch decom- 
position of the spectral problem. To this end, we follow [Gj [S] and define the one-parameter 
family of linear operators, referred to as Bloch operators, by 

L^u] := e-* x L[u]e* x , £e [-7r,vr], 

operating on L 2 er ([0, 1]). The L 2 spectrum of the linearized operator L[u] is readily seen to 
be given by the union of the spectra of the Bloch operators. By continuity of the spectrum 
on the Floquet parameter £, and the discreteness of the spectrum of the elliptic operator 
L[u] on the compact domain [0,1], it follows that the spectra of L[u] may be described as 
the union of countably many continuous surfaces //(£)• 

Continuing with this functional setup, we recall that any function g E L 2 (M) admits an 
inverse Bloch-Fourier representation 

9{x) = h f e * x 9& x W 

where x) = J2jez e 2lT ^ x g(S, + 2irj) is a 1-periodic functions of x, and g(-) denoting with 
slight abuse of notation the usual Fourier transform of the function g in the spatial variable 
x. Indeed, using the Fourier transform we have 

/oo rn rn 

e*^(0de = E / e i( -* +2 * j)x 9(Z + WW = / e**g(£, x)d£, 
-OO j^z 1 ^ J —7V 

where the summation and integral can be interchanged for Schwarz functions g. By similar 
computations, it is also seen by Parseval's identity that the Bloch-Fourier transform g(x) — > 
g(£,x) is an isometry of L 2 (W), i.e. 

/7r rl 
/ \g(£,x)\ 2 dx d£ =: \\g\\ L ^{i^(x))- 
-7T J o 

Moreover, the Bloch-Fourier transform diagonalizes the periodic-coefficient operator L[u], 
yielding the inverse Bloch-Fourier transform representation 

(2.5) e L ^g(x) = i- £ e^e^g^ x)d£ 

relating the behavior of the linearized system to that of the diagonal operators ^[i*]- 

We now discuss the spectral stability of the underlying solution it (a;) in more detail. To 
begin, notice that by the translation invariance of (jl.ip the function u'(x) is a 1-periodic 
solution of the differential equation L[n]f = 0. Hence, it follows that \i = is an eigenvalue 
of the Bloch operator LoM- From the structure of the equation, it is natural to assume 
that translation generates the only null direction of the operator LoM- Indeed, notice 
by (H2) and the secular dependence of the periodic on the wave-speed c, that variation 
in the translation direction xo is the only generator of 1-periodic null-directions of the 
linearized operator given by Noether's theorem. With this in mind, following [S], we make 
the following natural spectral stability assumptions: 
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(Dl) jj, = is a simple eigenvalue of i>o[it]. (Recall that £ = corresponds to co-periodic 
perturbations.) 

(D2) Ra(L(:[u)) < -6\£\ 2 for some constant 9 > 0. 

Assumptions (D1)-(D2) correspond to "dissipativity" of the large-time bahavior of the lin- 
earized system and are often referred to as strong or diffusive spectral stability assumptions 



Remark 2.1. By standard spectral perturbation theory [K] , (Dl) implies that the eigen- 
value //(£) bifurcating from /i = at £ = is analytic at £ = 0, with ^(£) = + + 
0(|£| 3 ), from which we find from the necessary stability condition 5fy/(£) < that K/Ui = 
and 3fyi2 < 0. Assumption (D2) thus amounts to the nondegeneracy condition 3fyt2 ^ 
together with the strict stability condition fficrL^ < for £ ^ 0. Condition (D2) is never 
satisfied in the scalar case n = 1, by Sturm-Liouville considerations, hence the emphasis in 
the title on systems of reaction diffusion equations. 

The goal of our analysis is to prove that the above spectral stability assumptions imply 
nonlinear L 1 f] H — > H K stability of the underlying periodic traveling wave. To this 
end, we use this spectral information to obtain bounds on the linearized solution operator 
e Z[«.]t_ ^ s we w jrj see ^ assumptions (D1)-(D2), along with (H1)-(H2), imply that the solution 
operator decays (in a suitable sense) polynomially in time at a fast enough rate to prove 
the nonlinear stability of the underlying spatially periodic solution. These bounds are 
established in the next section, after which we present our nonlinear iteration scheme. 

3 Linear Estimates 

In this section, we make use of the spectral stability assumptions of the previous section in 
order to prove bounds on the solution operator S(t) := e^"'*. These linearized estimates 
form the crux of the nonlinear analysis, presented in the next section. To begin, notice 
that by standard spectral perturbation theory [Kj, assumption (Dl) implies that the total 
eigenprojection -P(£) onto the eigenspace of L^[u] associated with the eigenvalue ^(£) bi- 
furcating from the (£,/•*(£)) = (^>0) state is well defined and analytic in £ for £ sufficiently 
small, since the discreteness of the spectrum of L^[u] implies that the eigenvalue is 
separated at £ = from the remainder of the spectrum of Lq[u]. In particular, there exists 
an eigenfunction (/(#,£) bifurcating from the q(x,0) = u'(x) state defined for [£| 1 such 
that 



[JZllOZllS]. 



(3.1) 



Li[u]q(x,£) = v{£)q{x,Q 



where, by assumption (D2), the function //(£) satisfies the estimate 



(3.2) 



»(m(0) < -m 2 



for some constant 6 > 0. 
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Our strategy is to treat the high- and low-frequency parts of the full solution operator 
S(t) separately since, as is typical, the low-frequency analysis is considerably more delicate 
than the corresponding high-frequency analysis. To this end, we introduce a smooth cut off 
function such that 

1, if |£| < e 
0, if Id > 2e, 

where e > is a sufficiently small parameter, and we split the solution operator S(t) into 
its low-frequency part 



HO 



(3.3) S^tMx) := -L J* e^mme^m, 

and high frequency part 



S H (t)g(x) :-- 



1 

2^ 



(1- 0(0^(0) e LfMt S(£,*K, 



by which one may readily check that S(t)g = (S 1 (t) + S H (t)) g by ([S3]) . As the low- 
frequency analysis is more delicate, we begin by deriving I? — > L p bounds on S n (t). 

Using the fact that L^[u] is a sectorial operator, and the spectral separation of from 
the remaining spectrum of £^[tt], standard semigroup theoryjHeJ [Pa] trivially implies that 
the bounds 



(3.4) 



e L S [u]t {l _ m p {0)dxg 

dxe L,W {l _ m p {0)g 



L2([0,l]) 



~ e ll5lk 2 ([o,i])j 



L 2 ([0,1]) 



5lli a ([o,i])j 



L 2 ([0,1]) 



e !I5||l 2 ([o,i])> 



for all t > and some constant 6 > 0. Together with (|2.4p . this yields immediately the 
following estimate. 

Proposition 3.1. Under assumptions (H1)-(H2) and (D1)-(D2), there exists a constant 
9 > such that for all 2 < p < oo and t > we have 

ll^(i)ffll^( H )<i^ (1/2 - 1/p) e- ei |b||^( R) - 
Proof. First, notice the bounds in (|3.4p and the triangle inequality imply that 



d™S n {t)g\\ LH 



< 



d™(i-mno)e L([u]t m,-) 



L 2 (x;[0,l}) 



< t -m/2 e -et 



ls(f> ■)IIl 2 (x-;[0,1])^ 



+-m/2-8t n II 
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for either m = or m = 1, where the final equality is justified by (|2.4p . thus justifying the 
first claim by taking m = 0. 

To prove the second inequality, note when p = oo Sobolev embedding and the above 
L 2 I? bound implies that 

\\S U $)9\\v»<n < {\\S n (t)g\\LHm) ■ \\d x S n (t)g\\ LHxm ) 1/2 < t-^e^Wgh^. 

The result for general 2 < p < oo now follows by LP interpolation. □ 

In order to analyze the low-frequency part of the solution operator S(t), we find it 
convenient to introduce the Green kernel 

G 7 (x,t;y) := S^t^x) 

associated with S 1 , and 

[Gl(x,t;y)] := ^)P^)^ [Sy(x)] 

the corresponding kernel appearing within the Bloch- Fourier representation of G 1 , where 
the brackets [•] denotes the periodic extensions of the given function onto the whole real 
line. Our first goal is to provide a useful representation for G 1 which incorporates the 
spectral properties (D1)-(D2) of the previous section. To begin, we introduce/recall a bit of 
notation. We denote by q(x,£) and q(x,£) the right and left eigenfunctions of the operator 
L^[u], respectively, associated with the eigenvalue bifurcating from the (£, /z(£)) = (0, 0) 
state. Moreover, we assume the normalization condition (q(-, £), q(-, S,))l 2 ([o 11) = !■ Then 
we have the following representation for the Green kernel G 1 . 

Lemma 3.2. Under assumptions (H1)-(H2) and (D1)-(D2), we have 

[G{{x,t;y)] =m^ )t q(x,Qq{y,Zl 
(3.5) G\x,t- iy ) = J- / e^-v) [G|(a;,t;y)] d£ 

Z7r JR 

= J_ f e*^U(.Z)e!*&q(x,Z)q(y,OdZ. 

Proof. The first equality is immediate from the spectral decomposition of elliptic operators 
on compact spatial domains. Moreover, using the fact that the Fourier transform (either 
continuous or discrete) of the delta function is unity we have 

= e-Xv^dteifc-v) = e -*v[6 y (x)]. 
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It follows that 



G I (x,t;y) 



2tt 



2tt 



which yields the second equality by recalling that cj) is supported in [— it, it]. 



□ 



In order to obtain nonlinear stability in the present context, it turns out that one can 
not simply use bounds as in the previous section on the function G 1 : see Remark 13. II below. 
Instead, we must take extra care in obtaining our low-frequency linearized estimates by 
separating out the slow-decaying translation mode from the faster-decaying "good" part of 
the solution operator. To this end, we define the function 



e(x,t;y) := -L / e*W<j>@e!*&q(y,OdZ 



and notice that 

G^x^y) = u'(x)e{x,t;y) + -L f e *(*-v) 0(^)^(0* ( q ( x ,£) - u'{x)) q(y,0^. 



2vr „ 

By the analyticity of q(x,£) in the variable £, we have that q(x,£) — u'(x) = 0(|£|) for 
|£| <C 1 and hence we expect the difference G I (x,t;y) — u'(x)e(x,t;y) to decay faster than 
the full Green kernel G . This is the content of the following proposition. 

Proposition 3.3. Under the hypothesis (H1)-(H2) and (D1)-(D2), the low-frequency Green 
kernel G 1 can be decomposed as 

G J {x, t; y) = u'(x)e(x, t; y) + G T (x, t; y) 

where, for all t > 0, and 2 < p < oo we have the estimate 



(3.6) 



sup 

y 



G'i-^y) <(1 + *) - 

Lp(R) 



§(i-V*>)-4 



Moreover, for all t > 0, 0<j,k,r,j + r<K + l we have the bound 



(3.7) 



sup 

y 



did«d r y e(x,t;y) 



< 



LP(x;. 



(! + *)• 



i(i_i /p) _(i+M 



Proof. Using the analyticity of the function q(x, £) on the variable £, we have that 
G\x,t-y) = ±- [ e^-^me^Oimiv.m 
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and hence the triangle inequality and (D2) yields 



sup 

y 



<?'(-, t;i/) 



< 



Moreover, noting that (|3.5p may be viewed itself as a Bloch-Fourier decomposition with 
respect to the variable z := x — y, with y appearing as a parameter, we may use (|2.4h to 
estimate 



sup 

y 



i$ sup 

V 



£ sup 

V 



me 
<(i+*r 3/4 



0(Oe^(g(x,e)-n'(x))g(y,e) 

lei 



L2(x;L2(0) 



'l«| 2 *| 



sup||g(y,-)||L°°( 5; h^]) 



where we have used in a crucial way the boundedness of g. By LP interpolation then, we 
obtain the desired LP bounds on G (x,t;y). Similar calculations yields the correspond- 
ing bound on e(x,t;y) by noting that y-derivatives do not improve decay while x- and 
t-derivatives improve decay by a factor of t~ l l 2 as above. □ 

Remark 3.1. Lt is important to note that the Green kernel G 1 does not decay fast enough 
in this one- dimensional setting to close our nonlinear iteration argument presented in the 
next section. Lndeed, using calculations as above (see also one can verify that 



sup ||G 7 (-,t;?/) 



\LP( 



for all t > and 2 < p < oo. Thus, by factoring out the translation mode from G 1 we 
gain an extra t -1 / 2 decay, which will end up being sufficient to close our nonlinear iteration 
arguments. 

We now combine the various bounds derived above to obtain estimates on the full Green 
kernel G(x, t; y). First off, let \{t) be a smooth cut off function defined for t > such that 
x (t) = for < t < 1 and x(t) = 1 for t > 2 and define 

e(x,t;y) := x(t)e(x,t;y). 

Corollary 3.4. Under the hypothesis (H1)-(H2) and (D1)-(D2), the Green kernel G(x,t;y) 
decomposes as 

G(x, t; y) = u'(x)e(x, t; y) + G(x, t; y) 
where for all t > 0, 2 < p < oo, and < j,k,r,j + r < K + 1, we have 



(3.8) 



G(;t;y)g(y)dy 



< (1+t) -i t 4(l/2-l/ P )| 



Li>( 
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and 
(3.9) 



did*d r v e(;t;y)g(y)dy 



< 



(1+t 



j+k 
2 



Proof. The bound (13, 8p follows immediately by considering the cases < t < 1 and t > 1 
separately. Indeed, using the L 2 — > L p high-frequency bounds in Proposition 13.11 for short 
time and the L p bound (|3.6p of Proposition 13.31 together with the triangle inequality for 
large time yields the desired result. The bound (|3.9p follows similarly, using the LP bound 
(13. 71) of Proposition 13.31 together with the triangle inequality. □ 



4 Nonlinear Stability 

With the above linearized estimates in hand, we are in a suitable position to prove nonlinear 
stability of the periodic traveling wave u{x) of the system of reaction diffusion equations 
(jl.ip . Our main result is as follows. 

Theorem 4.1. let u be a periodic standing-wave solution of (II. ip and let u(x,t) be any 
solution of (jl.ip such that \\u — u\\ L i^ nH K^ is sufficiently small. Then assuming (Hl)- 
(H2) and (D1)-(D2), there exists a constant C > and a function if)(-,t) € W K '°°(M.) such 
that for all t > and p > 2 we have the estimates 

\\u - u{- - i>(;t))\\ LP{R) (t) < C(l + tyfr-Vri-zWu^t) - u\\ LlnHK{R) \t =0 , 

3 

(4.1) - u(- - if}(-,t))\\ H K (R) (t) < C(l + tyi\\u(-,t)-u\\ L i {R)nH K {R) \ t =o, 

3 

\\(ipt,ipx)(-,t)\\ H K (R) < C{1 + ty±\\u(-,t) - u\\ L ^(K)nH K (M)\t=o, 

and 

(4.2) \\u-u\\ LP{R) {t), U{-,t)\\ LP(R) < C{l + ty^ {1 - 1/p) \\u(-,t) -u\\ L i (R)nH K {R) \t= . 
In particular, u is nonlinearly asymptotically L 1 n H K — > H stable with estimate 

(4.3) \\u-u\\ H K (R) {t), \\ip{-,t)\\ H K (R) <C(l + t)~^\\u{-,t)-u\\ L i mnH K iR) \ t=0 
for all t > 0. 

Remark 4.1. The regularity requirements stated here can be reduced to / £ C 2 and 
smallness of the initial data vq := u — u in L l n L°° as described in Remark 14.21 This is to 
be compared with Schneider's [S] assumptions / G C 4 and Vq small in a weighted H l l 2+5 
space (5 > 0) bounding ||(1 + |x| 2 )i;o||l°°, hence IbollLinL 00 ) by Sobelev embedding. 
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4.1 Nonlinear Perturbation Equations 

Let u(x, t) be a solution of the system of reaction diffusion equations 

u t = u xx + f(u) + cu x 

and define u(x, t) = u(x+ip(x, t),t) for some unknown function ip : R 2 — >■ R to be determined 
later. Moreover, let u(x) be a stationary solution and define 

(4.4) v(x, t) = u(x, t) — u(x) = u(x + ip(x, t),t) — u(x) 

Lemma 4.2. For v, u as above, we have 

u t - u xx - f(u) - cu x = (d t - L) u'(x)ip(x, t) + d x R 

+ (dt - d 2 x )S + (f'(v(x, t) - u(x)) - df(u(x))) ^, 

where 



r\ ( \ i/i / I , I / w , ( \ U x\ + K 



R := vipt -vip xx + (u x (x) + v x (x,t)) - -j-^- = O y\v\(\ip t \ + IV'xxl) + y l - \ifj \ ) 
and 

S:=-viJ> x = 0(\v\-\il> x \). 

Proof. Using the fact that u t — u xx — f(u) — cu x = 0, it follows by a straightforward 
computation that 

(4.6) ut + f(u) - u xx - cu x = u x ip t - u t tp x - (u x ^ x ) x + f(u)ip x , 

where it is understood that the argument of the function u and its derivatives appearing 
on the righthand side are evaluated at (x + Moreover, by another direct calculation, 

using the fact that 

L(u'(x)) = (a 2 + cd x + df(u)) u'(x) = 
by translation invariance, we have 

(dt - L) u'(x)ip = u x ip t - (u x ip x ) x - (cu x + u xx )tp x = u x ip t ~ (u x ipx)x + df(u)ip x . 

Subtracting, and using the facts that, by differentiation of (u + v)(x,t) = u(x + ip,t), 

U x +V x = U x (l + 1p x ), 
U t + V t = U t + u x ip t , 

so that 

ipx 



u x - u x - v x = -(u x + v x 



(4.8) 1 + ^ 

u t -ut-v t = -(u x + v x ) , 

1 + Ipx 
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we obtain 

ut + f(u) - u xx = (d t - L)u'(x)ip + v x ip t - v t Tp x - (v x tp x ) x 

+ ((u x + v x )-^-] +(df(v + u) -df(u))ip x , 

yielding (gJ3) by v x ip t - v t tp x = {vi/j t ) x ~ {vip x )t and (v x tp x ) x = (vip x ) xx - (vip xx ) x . □ 
Corollary 4.3. The nonlinear residual v defined in ()4.4p satisfies 

(4.9) (d t - L) v = (d t - L) u'( Xl )iP + Q + R x + (d 2 x - d t )S + T, 
where 

(4.10) Q := f(v(x, t) + u(x)) - f(u(x)) - df(u(x))v = 0(\v\ 2 ), 

(4.11) R := vipt - vip xx + (u x + v x ) 



1 + ipx ' 



(4.12) S:=wl> x = 0(\v\\il> x \), 
and 

(4.13) T := (df(v + u)- df(u)) ^ = 0(\v\\^ x \). 
Proof. Straightforward Taylor expansion comparing (|4.5p and u t — f(u) 

^^XX 

= 0. □ 

4.2 Integral Representation/ V>Evolution Scheme 

Using Corollary 14.31 and applying Duhamel's principle we obtain the integral (implicit) 
representation 



v(x,t) = u'(x)ip(x,t) + J G(x,t;y)v (y)dy 

r-T poo 



DC 



Pi l-OO 

+ / / G(x, t-s;y)(Q + R y + (8 2 - 8 S )S + T) (y, s)dy ds 

JO J-oc 

for the nonlinear perturbation v. Thus, if we define tp implicitly via the formula 



poo 

t/j(x,t) := - / e(x,t;y)v (y)dy 

-oo 

ft fOO 



pt poo 

/ / e(x,t-s;y)(Q + R y + (d 2 -d s )S + T)(y,s)dyds, 

JO J-oo 
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we obtain the integral representation 



(4.14) 



fOO ____ 

v(x,t)=l G(x,t;y)v (y)dy 

-oo 

r-T f-oo 



f l poo 

+ / / G(x, t-s;y)(Q + R y + {d 2 y - d s )S + T) (y, s)dy ds. 

JO J -oo 



Moreover, differentiating and recalling that e(x, t;y) = for < t < 1 we obtain 
(4.15) 



roo 

d*d™rP(x,t) := - I d?d™e(x,t;y)v (y)dy 



OO 

/ / $d?e(x,t-s;y)(Q + R y + (a$-d s )S + T)(y,8)dyd8. 

JO J -oo 



'0 J -oo 

Together, these form a complete system in the variables (v,dt'ip,d™ip) , < k < 1, < 
m < K + 1. In particular, given a solution of the system we may afterward recover the shift 
function ijj. 

Now, from the original differential equation (|4.9p together with (|4.15p , we readily obtain 
short-time existence and continuity with respect to t of solution (v,ifjt,ip x ) £ H K by a 
standard contraction-mapping argument treating the linear df(u)v term of the lefthand 
side along with Q, R, S, T, iftu' terms of the righthand side as sources in the heat equation. 

4.3 Nonlinear Iteration 

Associated with the solution (u,^t,tp x ) of the integral system (|4.14p - (|4.15p . we define 
(4.16) V (t):= sup \\(v,^ t ,i> x )\\ H K {xm (s)(l + s) 3 /\ 

0<s<t 

By short time H (R) existence theory, the quantity \\(v, ipt,4'x)\\H K (R) ^ s continuous so long 
as it remains small. Thus, r/ is a continuous function of t as long is it remains small. We 
now use the linearized estimates of Section 3 to prove that if r] is initially small then it must 
remain so. 

Lemma 4.4. For all t > for which rj(t) is sufficiently small, we have the estimate 

r ] {t)<C{E + 11 {tf) 

for some constant C > 0, so long as Eq := \\v(-, 0)||£i(R)nir K (:R) ^ s a ^ so sufficiently small. 

Proof. To begin, notice that by the descriptions of Q, T, R, and S in Corollary 14.31 we have 
that 

\\(Q,R x ,T)(-,t)\\ LHR)nL 2 {R) < ||(w ) v a;J ^iV'*)llii( a; . R )(*) + \\(v,v x ,iJt,iPx)\\ 2 L 2 (x . M) (t) 

< v (tf(i + t)-^ 
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so long as \\(v x ,ip x )(-,t)\\Loo^ — \\( v i'4>x)\\H K (x;R)(t) — ? ?(*) remains bounded, and likewise 

\\(d t - d 2 x )S(;t)\\ LHR)nL2m < \\(v^ x )f HHxm (t) + \\(v^ x )\\ 2 H 2 {xm (t)<v(t) 2 (l+ty 3 / 2 . 

Thus, applying the bound (|3.8p of Corollary 13.41 to representations (|4.14f) - (|4.15j) . we obtain 
for any 2 < p < oo the bound 

IK,i)llL,(M)<(l + i)^ (1 - 1/pM/2 ^o 

+ + t _ s) -3/4 {t _ a) -i(i /2 -l/p) (1 + s) -3/2 ds 

(4.17) 7o 

< (£ Q + (1 + ^-min(I(l-l/ P )+I,I(l-l/ P )+l) 

<(E + V (t) 2 ) (l + ty^-m-h 

and similarly using (|3.9p we have 

||(^,^.)(-,*)ll^ + i„(R) < (l + t)-^ {1 - 1/p) i?0 

(4.18) +v(t) 2 [ (l+t-sy^M-^i + sy^ds 

Jo 

yielding in particular that \\(tpt, ' i Px)\\h k + 1 * s arbitrarily small if Eq and rj(t) areH thus 
verifying the hypothesis of Proposition 14.51 below. By the nonlinear damping estimate given 
in Proposition 14.51 therefore, the size of v in H (R) can be controlled by its size in L 2 (R) 
together with H K estimates on the derivatives of the phase function ip. In particular, we 
have for some positive constants 0\ and #2 

<e-^El + (E + V {tf) 2 {l + t)-^ 

<(£ + ^) 2 ) 2 (i + r 3/2 - 

This estimate together with (|4.18p in the case p = 2 completes the proof. □ 

Proposition 4.5. Assuming (HI) - (H2), let v(-,0) € H K (R) (for v as in ([PI) ) and 
suppose that for some T > the H K (Nl) norm of v and the H K+1 (M) norms of ipt(-,t) and 
ip x (-,t) remain bounded by a sufficiently small constant for all < t <T. Then there are 
constants 9\ and #2 such that 

IK.*)I&*(R) ~ e ^IK' )H^(R)+ /V^*-) (\\v(;s)\\l Hm + \\(ipt,^)(;s)\\ 2 HKm ) ds. 

J 

for allO<t<T. 

2 Note that we have gained a necessary one degree of regularity in if), the regularity of if) being limited 
only by the regularity of the coefficients of the underlying PDE (|2.1|) . 
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Proof. Subtracting from the equation (|4.6p for u the equation for u, we may write the 
nonlinear perturbation equation as 

(4.19) v t - df(u)v - v xx - cv x = Q + u x ip t - u t ip x - (u x ip x ) x + f(u)t/j x , 

where it is understood that derivatives of u appearing on the righthand side are evaluated 
at (x + ip(x, t),t). Using (14. 8p to replace u x and ut respectively by u x + v x — (u x + v x ) 1 ^ 

and ut + vt — (u x + v x ) q% , and moving the resulting vfip x term to the lefthand side of 
(|4.19p . we obtain 

(1 + ip x )v t - v xx - cv x = -df(u)v + Q + {u x + v x )ip t 

(4.20) _ / _ x 

+ \ {u x + V x ) X +f(u)lp x 
\ 1 + 1p x J x 

Taking the I? inner product in x of Ylj=o l+gf ~ a S ams t (|4.20[) . integrating by parts, and 
rearranging the resulting terms, we arrive at the inequality 

dt\H;t)\\ 2 HKm < -e\\d* +1 v(;t)\\h m + C (\H;t)\\ 2 H K m + \m^x)(;s)\\ 2 HKm ), 

for some 6 > 0, C > 0, so long as remains bounded, and \\v (-,t)\\ H K m\ and 

IK^tj ipx)(', t) \\h k + 1 (R) remain sufficiently small. Using the Sobolev interpolation Hflil^-jrm) < 
C* -1 ll^ +1 ^llz,2(] R ) + ^HsUi^m) for (7 > sufficiently large, we obtain 

9t\\v(-,t)\\ 2 HK (R) (i) < -0|K->*)IIjW(r) + C (lk(-,t)|li 2(IR) + IK^V'aOMIIiWfR)) 
from which the desired estimate follows by Gronwall's inequality. □ 



Proof of Theorem \4-l\ Recalling that rf(t) is continuous so long as it remains small, it 
follows by continuous induction using Lemma 14.41 that rj(t) < 2CEq for all t > provided 
that Eq < 1/4C and (as holds without loss of generality) C > 1, yielding by (|4.16p the 
result (|4.ip in the case p = 2. Similarly, using (|4.17|) and (|4.18|) . the result (|4.ip follows for 
any 2 < p < oo with uniform constant C > 0. 
Finally, notice that by (|4.4[) we have 

u(x, i) — u(x) = v(x — ip(x, t),t) + (u(x — ip(x, t)) — u(x)) 

and hence the size of u(x, t) — u(x) in L p or H is controlled by the corresponding size of 
the function (v + ifj)(x,t) in the respective norm. Therefore, using (|4.17p along with the 
estimate 

U(;t)h P (R) < E (l + t)-3d-i/p) + ^ [\i + t- S yha-Vp) (1 + a) -§ da 

J 



< 



(^ + Tj(t) 2 )(l + t)-*( 1 - 1 /p), 



which follows by (j3.9H for all 2 < p < oo, we obtain (14. 2p . Similarly we obtain (|4.3p . which 
completes the proof. 
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Remark 4.2. Integrating by parts in representation formulae (|4.14p - (|4.15p to exchange 
the d y and (dy — d s ) derivatives on R and S for — d y and [dy + d s ) derivatives on G, and 
noting that G y and (G yy + G S ) satisfy L p bounds equal to the bounds obtained on G times 
the time-integrable factor (1 +i) 1 / 2 i~ 1 / 2 , we could alternatively close the iteration entirely 
within LP spaces, avoiding the need for a nonlinear damping estimate and obtaining LP 
stability, p > 2 with the same rates given in Theorem 14.11 for data merely small in L 1 nL°°H 
However, this is special to the semilinear case, whereas the damping argument given here 
applies also to more general quasilinear equations. 
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